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Abstract — In this work, we propose a low-complexity variable 
forgetting factor (VFF) mechanism for recursive least square 
(RLS) algorithms in interference suppression applications. The 
proposed VFF mechanism employs an updated component re- 
lated to the time average of the error correlation to automatically 
adjust the forgetting factor in order to ensure fast convergence 
and good tracking of the interference and the channel. Con- 
vergence and tracking analyses are carried out and analyti- 
cal expressions for predicting the mean squared error of the 
proposed adaptation technique are obtained. Simulation results 
for a direct-sequence code-division multiple access (DS-CDMA) 
system are presented in nonstationary environments and show 
that the proposed VFF mechanism achieves superior performance 
to previously reported methods at a reduced complexity. 

Index Terms — Interference suppression, adaptive receivers, 
recursive least squares algorithm, variable forgetting factor 
mechanism, time-varying channels 



I. Introduction 

Wireless communication channels are dynamic by nature 
and present significant challenges to the design of receivers 
in, lEl- In particular, the estimation of the parameters of 
the channels and the receive filters of wireless receivers has 
received significant attention in the past years and is key 
to the performance of wireless systems |]3]- In such nonsta- 
tionary environments, specifically for the slow time-varying 
channels, adaptive techniques are of fundamental importance 
and encounter applications in many wireless transmission 
schemes such as multi-input multi-output (MIMO) H-Q, 
orthogonal-frequency division multiplexing (OFDM) [SI and 
spread spectrum systems ifTOl - lfTSl . 

The most common adaptive estimation techniques are the 
stochastic gradient (SG) and recursive least squares (RLS) 
algorithms 1211. The works in ||3|, ET\ have considered 
standard SG algorithms with fixed step sizes to implement the 
receiver The authors in ||221 . Il28l have proposed some variable 
step-size schemes to accelerate the convergence speed of the 
SG based adaptive filters. The RLS algorithm is considered 
as one of the fastest and most effective methods for adaptive 
implementation ||2T| . In nonstationary wireless environments 
in which users often enter and exit the system, it is impractical 
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to compute a predetermined value for the forgetting factor 
Therefore, the RLS algorithm needs to be modified in order to 
yield satisfactory performance in time-varying environments. 
In this regard, one promising technique that has been con- 
sidered is to employ variable forgetting factor (VFF) mecha- 
nisms to adjust the forgetting factor automatically lISTI - llZTl . 
The classic VFF mechanism is the gradient-based variable 
forgetting factor (GVFF) algorithm which has been proposed 
in [21]. The work in [211 . besides the recursive expressions 
to adapt the receive vector, uses another SG recursion to 
control the forgetting factor where the gradient with respect 
to the forgetting factor is obtained based on the instantaneous 
squared error cost function. In ll25l . the authors proposed a 
modified variable forgetting factor mechanism using the error 
criterion with noise variance weighting for frequency selective 
fading channel estimation. A VFF mechanism based on the 
Gauss Newton approach was proposed in ll24l to improve 
the fast time- varying parameter estimation. In ll26l . ll27l . a 
modified GVFF mechanism based on the gradient of the mean 
squared error (MSB) rather than on the gradient of the instan- 
taneous squared error was investigated. Another approach is to 
perturb the covariance matrix whenever the change is detected 
II23I . however, the algorithm becomes sensitive to disturbance 
and noise. The existing VFF mechanisms are mostly derived 
from the architecture of the classic GVFF algorithm iBTI . and 
the computational complexity of which is proportional to the 
receive filter length. 

In this work, we propose a low-complexity VFF mechanism 
for adaptive RLS algorithms applied to linear interference 
suppression in direct-sequence code-division multiple access 
(DS-CDMA) systems. The proposed VFF mechanism employs 
an updated component related to the time average of the error 
correlation to automatically adjust the forgetting factor in order 
to ensure good tracking of the interference and the channel. We 
refer to the proposed VFF scheme as correlated time-averaged 
variable forgetting factor (CTVFF) mechanism. Convergence 
and tracking analyses are carried out and analytical expres- 
sions for predicting the mean squared error of the proposed 
adaptation technique are obtained. It should be noted that the 
proposed scheme is general and can be applied to any wireless 
systems with the RLS algorithm. Therefore, the scheme can 
be very useful in parameter estimation applications to multi- 
antenna and OFDM systems. 

The paper is structured as follows. Section HI] briefly de- 
scribes the DS-CDMA system model and the design of lin- 
ear receivers. The adaptive RLS algorithm and the existing 
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gradient-based VFF mechanisms are reviewed in Section |III] 
The proposed CTVFF mechanism is described in Section |IV] 
Convergence and tracking analyses of the resulting algorithm 
and the analytical formulas to predict the steady-state MSE 
and the steady-state mean value of the proposed CTVFF 
mechanism are developed in Section |V] The simulation results 
are presented in Section |VT] Finally, Section IVIII draws the 
conclusions. 



II. DS-CDMA System Model and Design of Linear 
MMSE Receivers 



Detecting a desired signal in DS-CDMA systems requires 
processing the received signal in order to mitigate the inter- 
ference and the noise at the receiver The major source of 
interference in DS-CDMA systems is multiuser interference 
(MUI), which arises due to the fact that users communicate 
through the same physical channel with nonorthogonal signals. 
Multiuser detection has been proposed as a means to suppress 
MUI, increasing the capacity and the performance of CDMA 
systems ll9l- lfT8l . The linear minimum mean squared error 
(MMSE) receiver implemented with an adaptive algorithm is 
one of the most prominent techniques since it only requires 
the timing of the desired user and a training sequence. In this 
work, we focus on adaptive linear receivers as extensions to 
other receivers are straightforward. 

Let us consider the downlink of an uncoded synchronous 
binary phase-shift keying (BPSK) DS-CDMA system with 
K users, TV chips per symbol and Lp propagation paths. A 
synchronous model is assumed for simplicity since it captures 
most of the features of more realistic asynchronous models 
with small to moderate delay spreads. Let us assume that the 
signal has been demodulated at the mobile user, the channel 
is constant during each symbol and the receiver is perfectly 
synchronized with the main channel path. The received signal 
after filtering by a chip-pulse matched filter and sampled at 
chip rate yields an M-dimensional received vector at time i 



vector h(i) with the multipath components are described by 



K 

r{t) = {Akbk{i)Ckh{{) + rjk{{)) + n(z), (1) 

fc=i 



where M = N + Lp - 1, n(i) = [ni{i) ... nAf(i)]'^ 
is the complex Gaussian noise vector with zero mean and 
E[n{i)n^ {i)] ~ cr^I whose components are independent and 
identically distributed, where (.)^ and (.)^ denote transpose 
and Hermitian transpose, respectively, and E[.] stands for 
expected value. The user symbols are denoted by bk{i), where 
we assume that the symbols are independent and identically 
distributed random variables with equal probability from the 
set {±1}. The amplitude of user k is Ak, and the signature 
of user k is represented by — [ah{l) . . . ak{N)]'^ . The 
M X Lp constraint matrix C^- that contains one-chip shifted 
versions of the signature sequence for user k and the Lp x 1 



C, = 



afc(l) 

ak{N) 






afc(l) 



ak{N) 



hL^-i{i) 



(2) 

where the M x 1 vector Cfch(i) denotes the effective spreading 
code. The vector rif.(i) denotes the intersymbol interference 
(ISI) for user fc, here we express the ISI vector in a general 
form that is given by ?7j.(i) = Akhk{i - l)HPpfe + Akhkii + 
l)H'*pfc, where the M x N matrices H^* and account for 
the ISI from previous and subsequent symbols, respectively, 
and can be given as follows 



H 



p — 









ho{l + l) 



(^-l) ... hi{i-l) 


... 




(3) 



/iL,-2(« + l) ... haii + 1) ... 

The linear MMSE receiver design is equivalent to deter- 
mining a finite impulse response (FIR) filter Wfc(i) with M 
coefficients that provide an estimate of the desired symbol as 
follows 

- wf(i)r(i), (4) 

where the detected symbol is given by bk{i) = 
sign{¥t[wj^ {i)r{i)]}, where the operator 3?[.] retains the 
real part of the argument and sign{.} is the signum function. 

The design of the receive filter Wfc(z) is based on the 
optimization of the MSE cost function 



^MS£;(wfc(i)) = E \bk{i) 



{i)v{i)\ 



(5) 



By minimizing (|5]l, the MMSE receive filter is given by ||9] 



Wo 



R ^s, 



(6) 



where R = E[v{i)Y" {i)] = ^f/^-^ A2,C^yh(i)h^(i)C^ + 
Ef^=i4,H^P,.'P^H^'^ + Ef^.iA2,H^p,,p^ff^ + a2l 
and s = £'[6^(i)r(i)] = AfcCfch(i). The minimum value of 
the mean squared error is given by 



= 1 - A^h 



(^)CfR-^Cfch(^) 



(7) 



III. 



ADAPTIVE RLS ALGORITHMS AND PROBLEM 
STATEMENT 



In this section, we first describe the adaptive RLS algorithm 
to estimate the parameters of the linear MMSE receiver for 
multipath channels using a least-squares (LS) approach. Then, 
the existing gradient-based VFF mechanisms are introduced. 
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A. Adaptive Receiver 

Let us consider the time-averaged cost function 

i 

JLsil) - J2 ^'""l^feW - wf (z)r(n)|2 (8) 

n=l 

where A denotes the forgetting factor By taking the gradient 
of ([Sj with respect to w^(i) and setting it to zero, after 
further mathematical manipulations we have the adaptive RLS 
algorithm as follows II2TI 



where 



and 



Wfe(i) = Wkii - 1) + k(i)e*(i), 
R-i(i - l)r(i) 



A + r»(i)R-i(i--l)r(i) 
eW=6fcW-wf(i-l)r(z), 



(9) 

(10) 

(11) 



and the estimate R ^{i) is updated by 

IL-\i) = X-^Ti-\i - 1) - X-^k{i)r" {i)R-\i - 1). (12) 

The adaptive algorithm is implemented by using (l9]l- (fT2l i with 
appropriate initial values R^^(O) and Wfc(O). The algorithm 
has been devised to start its operation in the training (TR) 
mode, and then to switch to the decision-directed (DD) mode. 
The problem we are interested in solving is how to devise a 
cost-effective mechanism to adjust A, which is a key factor 
affecting the performance of RLS-based algorithms and the 
receivers. 

B. Gradient-based Mechanisms 

In order to adjust the forgetting factor automatically, the 
adaptive rule for the classic GVFF mechanism is derived by 
taking the gradient of the instantaneous cost function \bk{i) — 
wj^ {i)r{i)\'^ with respect to the variable forgetting factor A(i) 
||2T1 . The adaptive expressions are given by 



X{i)= A(i-l) + /i5R[ 



dX 



dX 



(l-k(z)r^(z)) 



dX 



dwkii-l) dR-^{i) 



(13) 



dX 



dX 



X-\^){I-^^)r"i^)) 



dX ' 



r(i)e*(i), 
(14) 

(15) 



where 



.•Ja- 



denotes the truncation to the limits of the range 



[A~, A+j, fi denotes a step-size. The adaptive receiver with the 
GVFF mechanism 1211 is implemented by using (l9])-(fT5ll with 
suitable initial values. By using the error criterion with noise 
variance weighting and the MSE rather than the instantaneous 
squared error, two modified GVFF mechanisms were recently 
proposed in 1251 and l27l . We refer to them as weighting 
GVFF (WGVFF) mechanism and mean squared error GVFF 
(MGVFF) mechanism, respectively. Fig.[T]illustrates the block 
diagram of the adaptive receiver with a variable forgetting 



factor mechanism. In the following section, we will describe 
the proposed low-complexity variable forgetting factor mech- 
anism. 



r(i) 



Receiver 



Adaptive 

RLS 
Algorithm 



VFF 
Mechanism 



Re(.) 



1' 



b(i) 



Fig. 1 . Structure of the adaptive receiver witli the variable forgetting factor 
mechanism 



IV. Proposed Variable Forgetting Factor 
Mechanism 

In this section, we first introduce the proposed low- 
complexity CTVFF mechanism that adjusts the forgetting 
factor of the RLS algorithm that equips the adaptive receiver. 
Then, we present the computational complexity analysis for 
the proposed CTVFF mechanism and the existing gradient- 
based mechanisms. 

A. Adaptive CTVFF Mechanism 

Based on the variable step-size mechanism for least mean 
square (LMS) algorithms in ||281 . the motivation for the 
proposed CTVFF mechanism is that for large estimation errors 
the adaptive algorithm will employ smaller forgetting factors 
whereas small estimation errors will result in an increase of the 
forgetting factor in order to yield smaller misadjustment. Thus, 
the CTVFF mechanism employs an exponential weighting 
parameter that controls the quality of the estimation error, and 
we have devised the following time-averaged expression 



7(^) = (5l7(^-l)+(52p'(^), 



(16) 



i)^h(erle-6(i^^)k^(I)afl(Jvl5^(z^HX^4»)4 5i is close to 1 and 
§2 is set to be a small value. The quantity p{i) denotes the 
time average estimate related to the correlation of — 1) — 
'w^ {i — l)r(i — 1)1 and \hk{i) — wj^ {i)r{i)\, which is given 
by 

= S.piz - 1) + (1 - 5s)\bkii - 1) - wf - l)r(z - l)\\bk{i) - 

(17) 

where < ($3 < 1, and 5^ controls the time average estimate. 
We set ($3 to be a quantity close to 1. 

The updated component p{i) is a small value, and it changes 
rapidly as the instantaneous value of the correlation between 

|5fe(i-l)-wf (i-l)r(i-l)| and - wf (i)r(i)|. The 
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use of 7(1) has the potential to provide a suitable indication 
of the evolution of the error correlation. Thus, unlike the 
conventional GVFF mechanism we aim to design a simpler 
mechanism that adjusts the forgetting factor automatically 
based on 7(1). Note that the forgetting factor should vary 
in an inversely proportional way to the value of the error 
correlation. We have studied a number of rules and the 
following expression is a result of several attempts to devise 
a simple and yet effective mechanism 

I A+ 



1 



(18) 



.l + 7(«). 

Indeed, the CTVFF mechanism is simple to implement and 
a detailed analysis of the algorithm is possible under a few 
assumptions commonly made in the literature. The proposed 
low-complexity CTVFF mechanism is given by ( fT6] l, ( fTTl ) and 
( fTSl l. The value of variable forgetting factor X{i) is close to 1, 
and it is controlled by the parameters Si, 62 and S3. We will 
investigate the sensitivity of the parameters in Section |VT] A 
large prediction error will cause the updated component 7(1) 
to increase, which simultaneously reduces the forgetting factor 
X{i) and provides a faster tracking. While a small prediction 
error will result in a decrease in the updated component 7(1), 
thereby the forgetting factor X{i) is increased to yield a smaller 
misadjustment. 

B. Computational Complexity 

We describe the computational complexity of the proposed 
CTVFF and the adaptive GVFF, WGVFF and MGVFF mech- 
anisms II2TI . ||25]| . II27I . Table H] shows the additional computa- 
tional complexity of the algorithms for multipath channels. We 
estimate the number of arithmetic operations by considering 
the number of complex additions and multiplications required 
by the mechanisms. The CTVFF and MGVFF mechanisms 
have constant computational complexity for each received 
symbol while the adaptive GVFF and WGVFF techniques 
have additional complexity proportional to the length M of 
the receive filter. An important advantage of the proposed 
adaptation rule is that it requires only ten operations. 

TABLE I 

Additional Computational complexity. 





Number of operations per symbol 


Algorithm 


Multiplications 


Additions 


CTVFF 


7 


3 


GVFF OD 


7M^ + 4M + 2 


7M^ + M 


WGVFF (25) 


7M^ + 4M + 9 


7M^ + M + 1 


MGVFF (27| 


29 


18 



V. ANALYSES OF THE PROPOSED ALGORITHMS 

In this section, we first show the convergence of the mean 
weight vector and derive the steady-state MSE expression for 
the adaptive receiver with the proposed CTVFF mechanism 
in the case of invariant channels. Then, we give the tracking 
analysis of the proposed scheme for the scenario of time- 
varying channels. Finally, the expression for the steady-state 
mean value of the CTVFF mechanism is derived. 



A. Convergence of the Mean Weight Vector 

We show the convergence of the mean weight vector of the 
proposed algorithm. Firstly, we give two expressions which 
will be used for the following derivation: 



and 



k(i) R-i(i)r(i) (19) 



lim R~^(i) w (1 - £;[A(oo)])R"^ (20) 



The proof of (fT9] l and ( |20] | are shown in the appendix. The 
relevant description is given in ll29l - ll3n . 
Let us recall ^ and assume 

- wfrW+Co(«), (21) 

where ^o(*) denotes the independent measurement error with 
zero mean and variance CTq, where CTq = 1 — AfcW^Cfch — 
Afeh^C^wo + WjJ'^Rwo, which is computed by ( I2TI 1. By 
defining e(i) = Wfe(i) — wq and using (|9]l we have 

e{i) = Wfc(i) - Wo = Wfc(i - 1) - + k(i)e*(i), (22) 

based on (fT9] l we have 

e(z) = e(* - 1) + K-\i)r{t){bUt) - r^(*)wfe(z - 1)) 
= e(* - 1) + R-i(z)r(z)(eo*(*) " ^" i^H^ - 1)) 
= (I - Il-\^)r{^)r"{^))e{^ - 1) + ■R-\t)r{t)Co{t). 

(23) 

When i ^ cxD, by substituting (|20] | into (l23T l we obtain 

e{i) = (I - (1 - E[X{oo)])R-^r{i)r" {i))e{i - 1) + (l - E[X{oo)])n- 

(24) 

We note that 1 — £'[A(oo)] is a small value. According to 
the direct averaging method described in 11211 . the solution 
of the stochastic difference equation (l24l i. operating under the 
condition of a small parameter 1 — £'[A(oo)], is close to the 
solution of another linear stochastic difference equation that 
is obtained by replacing the matrix (1 — £'[A(oo)])r(i)r^(i) 
with its ensemble average (1 — £^[A(oo)])R,. We may write the 
new expression as follows 

e{i) « (I - (1 - £;[A(oo)])R-iR)e(i - 1) + (l - E[X{oo)])R~^r{i)^* 
= E[X{oo)]eii - 1) + (1 - E[Xioo)])R-^r{i)CUi)- 

(25) 

By taking the expectation of ( |25] ). when i — cx), we have 

E[e{i)] « E[X{oo)]E[e{i - 1)]. (26) 

Since < £^[A(oo)] < 1, the expected weight error converges 
to zero. 

B. Convergence of MSE 

Then, we show the convergence of MSE for the proposed 
algorithm and give an analytical expression to predict the 
steady-state MSE. 
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Based on ( |25] ) we have 



By redefining e{i) = Wk{i) — wo(i) and following the 
aforementioned approach, we have 



E^[X{oo)]E[eii - l)e 



HfA^l-^^ e(i) = Wfe(i)-wo(i) 

, _ „ _ ^ 'vvhii — I) — "Woii — 1) — q(i) + 'k(i)e* (i) 

+ i;[A(w)](l-£;[A(c5o)])£;[e(i-l)]i;[eoWr^WR (34) 

(27) 



Note that when i — !> oo, i?[e(i — 1)] w 0, we have 

« £;2[A(oo)]e(i - 1) + (1 - E[X{oo)]YR-^al (28) 
Since < E^[X{oo)] < 1, we obtain 

■1 - E[X{oo)] 



Using ( |20] | when ?' becomes large, we obtain 

e{i) « (I - (1 - £;[A(oo)])R-ir(i)r^(i))e(i - 1) + (l - £;[A(oo)])R-^ 

(35) 



0(oo) 



R 



,l + £:[A(cx))] 
Note that the steady-state MSE is given by 

lim ^{i) = lim E[\bk{i) - wf (i - l)r(i)p] 

= lim (S(t) + 1 - ylfcW^Cfch - Afeh^Cf wo), 

(30) 

and 

S(z) = - 1) + Wo«)r(z)r^(z)(e(z - 1) + Wp)] 



Since 1 — E[X{oo)] is a small value, we may solve ( 135] ) for 
the weight error vector e{i) by employing the direct averaging 
method fZl\, which was mentioned previously. We obtain the 
(29) following expression 



e{i) « i?[A(^)]e(i - 1) + (1 - i?[A(oo)])R-ir(z)e*(i) - q(i). 

(36) 

Note that the vector q(i) is an independent zero mean vector, 
we have 

&i^) = E[ei^)e"i^)] 

« E^[X{oo)]@{i - 1) + (1 - S[A(oo)])'r-V2 + E[qii)q"{i)], 

(37) 

and when i becomes large, the MSE in a time-varying envi- 



w^Rwo + tr[R0(z - 1)] + S[e«(i - l)]i;[r(i)r^(i)W"^n" is given by 



+ E[w^r{i)r"{i)]E[e{t-l)]. 



(31) 



l)Rwo(j - 1) + tr[R0(i)] + 1 - Akw^{i - l)Cfch - / 



When i — > cx), we have ^(oo) = v\^,^Rwo + tr[R0(oo)]. 
Multiplying ( |29l ) by R we have 



(38) 



By using ( |37] |. when i becomes large, we have 



tr[R0(oo)] 



1 -£;[A(oo)] 

TT"e[amI 



(32) tr[R0(z)] 



1 -£;[A(oo)] 
l + £;[A(oo)] 



tr[Rg[q(»)q^(»)]] 

l-E^[X{oo)] 



(39) 

Therefore, the steady-state MSE in (jSO) can be approximated The MSE for a situation in which the adaptive receiver 
by the following expression 



1 - E[X{oo)] 



is tracking a channel can be computed by the following 
expression 



where 



1 + E[Xioo)] 
l-A2h^(z)CfR-iC.h(z). 



(33) 



e(00) « 



1 - -B[A(oo)] 
l + £;[A(oo)] 



tr[Rg[q(z)qg(z)]] 



(40) 

Note that we need to compute the quantities £'[A(oo)] and 
E[q{i)q^ (i)] to calculate the above expression. 
C. Tracking Analysis 

In a nonstationary scenario, the optimum solution takes D. Steady-State Mean Value of the CTVFF Mechanism 
on a time-varying form. It brings a new task to track the In order to derive the expression for the steady-state mean 
minimum point of the error-performance surface, which is no value of the CTVFF mechanism, we first show the conver- 
longer fixed. We investigate the convergence properties of the gence and derive the expression for the steady-state statistical 
proposed algorithm in this case. In time-varying channels, the property of the updated component 7(1). 
optimum filter coefficients are considered to vary according By recalling ([TtIi, the estimate of p{i) can be alternatively 
to the model wo(i) = wo(i — 1) + q(i), where q(i) denotes written as 
a random perturbation lISTI . We assume that q(z) is an in- 

dependently generated sequence with zero mean and positive ^(i) = (1 _ ^ _ ; _ 1) _ wf - I - l)r(i - I - \)\\bk{i 

definite autocorrelation matrix E[q{i)c^ {%)]. This is typical in 



the context of tracking analyses of adaptive filters 



(41) 
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and by squaring flTb we have 



i~l i-1 
n2 \ ^ \ ^ 



where is computed according to the Jakes' model ||2l. 

For the initial values, we set R~^(0) = I, Wfc(O) = 0.01 x 1, 



1=0 j=0 



/-l)-wf(*-Z 



(0) 



vector. The parameters of the MGVFF and WGVFF mecha- 

X \bk{i -3-1)- wf - J - l)r(^ - j - l)||6fc(^ - lfi^mm^mmM)f^- The simulations are averaged over 

1000(1 runs. We set the power of the desired user \Ai\ = 1. 



and 



I, where 1 denotes an all-one 



X \bk{i - j) - [i - j)r{i - j)\ 



(42) 



When i — > oo, we assume that \bk{i—l)—wj^{i~l)r{i — l)\ and 
\bk {i — j) — (i ~ j)r{i — j)\ are uncorrected, thus, we have 
E[M^-l)-^^^i^-l)vi^-l)\M^-J)-^^j^i^-J)r{^-J)\] = 

E[Mi-l)-^j^{i-l)r{t-l)\]E[Mi-j)-^ni-jMi- 
Ki 0, where Vj ^ I. Hence, the expectation of p^{i) is 

given by 



I 



Note that < 6-i 



rH 



< 1, by using 
we obtain 



lim, 



E 



E[p^{^)] 



2 

min 



(43) 

\bk{^) - 

(44) 



Note that < (5i < 1, based on (fTSl l we obtain the steady-state 
mean for 7(1), 



A. Effects of 5i, 82 and S3 

Fig.[2]shows the effect of 53 on the proposed CTVFF mech- 
anism. We show the received steady-state signal to interference 
plus noise ratio (SINR) versus S3 for different values of Si and 
52 under different scenarios. For the first scenario, we consider 
that the system has six users including two users operating at a 
power level 3dB above and one user operating at a power level 
6dB above the desired user's power level. The channel fading 
l)r{i^li^i^f]B[\lk^i^l^-J^^^W^fh{^'^P}\%rome is given 
by po = OdB, pi = —6dB and p2 = —lOdB. For the second 
scenario, the system has ten users which have the same power 
level. In this case, we consider a static channel, the channel 
parameters of which are given by Hq = OdB, hi = -~3dB and 
/i2 = —GdB. In Fig. |2] we can see that the performance does 
not change too much as the value of 63 varies for different 
values of Si and S2. We set 7- = 0.98 and 7+ = 0.99998, 
and use 15dB for the input signal to noise ratio (SNR). We 
tuned S3 — 0.99 for the following simulations. 



C2 



(45) 



(l-Si)il + S3)' 

From ( fTSI l and ([TtI i. we can see that the quantity 7(1) is a 
small value, and varies slowly around its mean value. 
Thus, when i — > cx), by using ( fTSl ) and (l45l l we assume that 
the steady-state mean value of the variable forgetting factor 
for the CTVFF mechanism is given as 

1 



E[Xi^)] 



1 + E[j{oo)] 

{1-Si){l + S3) 



(46) 



il-Si)il + S3) + S2il-S3)e. 

The above value can be substituted in ( [33] l and (|40l i to predict 



the MSE. 



VI. Simulations 



In this section, we first adopt a simulation approach and 
conduct several experiments in order to verify the effectiveness 
of the proposed CTVFF mechanism in RLS algorithms applied 
to interference suppression problems with DS-CDMA systems. 
Then, we verify the effectiveness of the proposed analytical 
expressions in ( [33] l, (|40l i and ( |46] l to predict the performance 
of the adaptive linear receivers in several situations. The DS- 
CDMA system employs Gold sequences as the spreading 
codes, and the spreading gain is = 15. The sequence 
of channel coefficients for each path is given by hf{i) = 
Pfaf{i){f = 0,1,2). All channels are normalized so that 



E 



1, 



(47) 



N=15, SNR=15clB 



ii q g B U ^ M S» ^ < i ^ 



fading channel f T=0.00001 , K=6 



static channel, K^10 



-e-s 


=0.96, 82=0.005, K=6 




=0.9879, 82=0.001, K=6 




=0.9819, 8^=0.001, K=10, static 




=0.9819, 8^=0.0005, K=10, static 




=0.93, 82=0.0005, K=10, static 




=0.9759, 82=0.001, K=6 



0.95 
8, 



Fig. 2. 

SNR -- 



Steady-state SINR versus ^3 
15 dB. 



for different values of 5i and 82 



We investigate the effects of Si and S2 on the CTVFF 
mechanism, we show the received steady-state SINR versus 
Si for S2 = 0.015, 0.005, 0.001, 0.0005, 0.0001 under different 
scenarios, where we tuned S3 ~ 0.99. We use 15dB for 
the input SNR, and set 7- = 0.98 and 7+ = 0.99998 for 
the CTVFF mechanism to guarantee the stability. The results 
shown in Fig. [3] (a) and (b) are based on a static state channel. 
The time-invariant channel parameters are given by hg = OdB, 
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hi = —3dB and h2 = —6dB. In Fig. [3] (a), the system 
has six users with the same power level. In Fig. |3] (b), the 
system has ten users with the same power level. It is observed 
that, firstly, the optimum choice of the pair (61,62) is not 
unique, for each value of 62 we can find a relevant 61 to have 
the best performance. Secondly, with the increasing of 62 the 
performance becomes less sensitive to 61. Moreover, when 62 
increases to 0.005, the performance becomes insensitive to 
61. From the results, we can obtain that 61 = 0.9819 and 
62 = 0.005 are one of the best choices for the scenario in Fig. 
|3] (a) and Si = 0.9879 and 62 = 0.001 are one of the best 
choices for both scenarios, which shows that it is possible for 
us to find a pair {61, 62) that works well for different scenarios. 



specific operating point, and stored at the receiver. 



(a) N=15, K=6, S =0.99, f T=0.00001 



(b) N=15, K=10, 5 =0.99, f 7=0.00001 



(a) N=15, K=6, 6 =0.99 



(b) N=15, K=10 ,8 =0.99 





and 



index 



Steady-state SINR versus 5i for different values of 52. (a) K = 6, 
3dB and one 6dB high power level interferers, (b) K = 10, with four 
two 6dB high power level interferers. f^T = 0.00001, S3 = 0.99. 



Fig. 3. Steady-state SINR versus 5i for different values of 52. (a) K = 
6, with equal power level interferers, (b) K = 10, with equal power level 
interferers. 53 = 0.99, static channel. 

The results in Fig.|4](a) and (b) show the steady-state SINR 
versus 6i for different values of ^2 based on the channel with 
fdT = 0.00001. The channel has a power profile given by 
Po = MB, pi = -MB and ps = -IQdB. We tuned 5s. = 
0.99. In Fig. H] (a), the system includes six users including 
two users operating at a power level 'idB above and one user 
operating at a power level QdB above the desired user's power 
level. In Fig. |4] (b), the system has ten users including four 
users operating at a power level idB above and two users 
operating at a power level QdB above the desired user's power 
level. We can see that the previous findings on effects of 6i 
and 62 still hold for the experiments. In particular, in Fig. 
|4] (a), when 62 increases to 0.015, the performance becomes 
insensitive to 61. In this case, we found that 61 = 0.934 and 
62 = 0.005 are one of the best choices for both scenarios. 
In the following experiments, for a given scenario we first 
choose the optimum parameters, and then fix them through 
the simulations. In practice, these optimized values should be 
obtained by experimentation for a given range of values or a 



o.offlj \ Pcrfonnance of the Proposed Algorithm 

We: first choose the received SINR as the performance 
to evaluate the convergence performance in nonsta- 
tionat-y scenarios. In this simulation, we assess the SINR 
^perfor mance of the proposed CTVFF mechanism, the adaptive 
GVFF mechanism, the MGVFF mechanism, the WGVFF 
mechanism, the fixed forgetting factor mechanism and the 
adaptive SG receiver The results shown in Fig. |5]illustrate that 
the performance in terms of SINR of the analyzed algorithms 
in a nonstationary scenario. The channel has a power profile 
given by po = MB, pi = —MB and p2 — —\MB. The 
system starts with six users including two users operating at a 
power level 3 dB above and one user operating at a power level 
6 dB above the desired user's power level. At 1000 symbols, 
four interferers including two users operating at a power level 
3 dB above and one user operating at a power level 6 dB 
above the desired user's power level enter the system. The 
normalized Doppler frequency is fdT = 1 x 10~^. From Fig. 
|5] we can see that the proposed CTVFF mechanism with the 
adaptive receiver achieves the best performance, followed by 
the WGVFF mechanism with the RLS receiver and the GVFF 
mechanism with the RLS receiver The MGVFF mechanism 
does not work well in the nonstationary scenario of multipath 
fading channels. The algorithms process 250 symbols in the 
TR mode. For the simulation, we tuned the forgetting factor 
A = 0.997 for the fixed forgetting factor scheme. For the 
GVFF mechanism, we tuned A" = 0.992, A+ = 0.99998, 
A(0) = 0.998 and /i = 0.0025. For the proposed CTVFF 
mechanism, we tuned A" = 0.98, A+ = 0.99998, 61 = 0.934, 
62 = 0.005, 6i = 0.99, 7(0) = and p(0) = 0. 



8 



N=15, SNR=15dB, f T=0.00001 




has 
nd 



a power profile given by pq ~ OdB, pi = —idB 
P2 = —6dB. The normalized Doppler frequency is 
= 1 X 10^"*. The system starts with eight users including 
a power level 3dB above and one 



200 400 600 800 1000 1200 1400 1600 18 

number of received symbols 



Fig. 5. SINR performance in nonstationary environment of multipath fading 
channels. SNR= 15dB. f^T = 0.00001. 



The result shown in Fig. |6] illustrates the variation of the 
forgetting factor value of the proposed CTVFF mechanism 
versus number of received symbols in the experiment of Fig. 
|5j We can see that the forgetting factor value becomes small 
at the first several iterations due to the large prediction error 
and it increases gradually until the system converges. At 
1000 symbols, when the new interferers enter the system, the 
forgetting factor will be recomputed. 



-A- RLS fixed 

* RLS GVFF 

□ RLSWGVFFr. 

V RLS CTVFF three users operatmg at 

O MMSE-SG [iser operating at a power level 6dB above the desired user's 

X RLSMGVFF \ , ° , , .c • ^ • , j- 

— — —power level. At 1000 symbols, tour mterterers mcludmg two 

users operating at a power level 3dB above the desired user's 
power level enter the system. From the results, we can see 
that the adaptive receiver with the CTVFF mechanism works 

^11 ^ 611 and still outperforms the conventional techniques. For the 

proposed CTVFF mechanism, we maintained the parameters 
n the simulation of Fig. |5] We tuned the forgetting factor 
A = 0.995 for the fixed forgetting factor scheme. For the 
(jVFF mechanism, we tuned A" = 0.992, A+ = 0.99998 and 
/t = 0.0035. 

N=15, SNR=15dB, f T=0.0001 



^^^^^^^^^ 






— — ex _ 


^ — 



N=15, SNR=15dB 



1.001 



0.999 



3 0.998 - 



a 0.997 - 



pi 0.996 



0.995 



0.994 - 



0.993 




1000 

number of received symbols 



Fig. 6. Forgetting factor variation of the proposed CTVFF scheme in 
nonstadonaiy scenarios. SNR= 15dB, /^jT = 1 X 10~^. 

In Fig. |7] we depict the SINR performance of the analyzed 
algorithms in a different nonstationary scenario. The channel 



600 800 1000 1200 1400 

number of received symbols 



7. SINR perfonnance in nonstationaiy environment of multipath fading 
chttnnels. SNR= 15dB. f^T = 0.0001. 



We examine the SINR performance of the CTVFF mech- 
anism in a static environment. We compare the CTVFF 
m;chanism with the adaptive GVFF mechanism, the MGVFF 
'schanism, the WGVFF mechanism, the fixed forgetting fac- 
tor mechanism and the adaptive SG receiver Fig. |8] illustrates 
th; SINR performance versus the number of received symbols 
for a static environment, where the time-invariant channel 
parameters are given by ho = OdB, hi = —3dB and /i2 = 
MB. The system includes six users which have equal power 
level. The results show that the proposed CTVFF mechanism 
works well in a static environment, and it outperforms the 
§|ggr conventional schemes. For the simulation, we tuned the 
forgetting factor A = 0.9995 for the fixed forgetting factor 
scheme. For the GVFF mechanism, we tuned A^ = 0.992, 
A+ = 0.99998, A(0) = 0.998 and ^ = 0.006. For the proposed 
CTVFF mechanism, we tuned A" = 0.98, A+ = 0.99998, 
Si = 0.9879, 62 = 0.001 and ^3 = 0.99. 

We show the bit error rate (BER) performance of the 
following algorithms as the fading rate of the channel varies. 
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N=15, K=6, SNR=15dB 



N=15, K=5, SNR=15dB 






RLS fixed 


* 


RLS GVFF 


□ 


RLS WGVFF 


V 


RLS CTVFF 


o 


MMSE-SG 


X 


RLS MGVFF 



500 1000 

number of received symbols 




RLS fixed \ 
-H- RLS GVFF 
-0-RLS WGVFF 
—I— RLS CTVFF 
-^SG 



Fig. 8. SINR performance in a static environment. SNR= 15dB. K - 



i.e., the RLS receiver with the proposed CTVFF mechanism, 
the RLS receiver with the GVFF mechanism, the RLS receiver 
with the WGVFF mechanism, the RLS receiver with the fixed 
forgetting factor mechanism and the adaptive SG receiver In 
this experiment, the values of the forgetting factors for the al- 
gorithms which employ fixed forgetting factors are optimized, 
and we tuned suitable forgetting factors for each value of 
fdT. The channel has a power profile given by po = OdB, 
pi = —6dB and p2 = —lOdB. In Fig. |9] we can see 
that, as the fading rate increases, the performance gets worse, 
and our proposed CTVFF scheme outperforms the existing 
schemes, which shows the ability of the RLS receiver with the 
CTVFF mechanism to deal with channel uncertainties. For the 
GVFF mechanism, we tuned A" = 0.993, A+ = 0.99998 and 
yU = 0.003. For the proposed CTVFF mechanism, we tuned 
A- = 0.98, A+ = 0.99998, Si = 0.988, S2 = 0.001 and 
63 = 0.99. 

Fig. [To] (a) and (b) illustrate the BER performance of the 
desired user versus SNR and number of users K, respectively, 
where we set fdT = 1 x lO^"*. We can see that the best perfor- 
mance is achieved by the adaptive receiver with the CTVFF 
mechanism, followed by the RLS receiver with the WGVFF 
mechanism, the RLS receiver with the GVFF mechanism, the 
RLS receiver with the fixed forgetting factor mechanism, the 
adaptive SG receiver and the conventional Rake receiver In 
particular, the adaptive receiver with the CTVFF mechanism 
can save up to over 4dB and support up to two more users in 
comparison with the receiver with the WGVFF mechanism, at 
the BER level of 10^^. 1500 symbols are transmitted and 250 
symbols are used in TR. 

C. Convergence and Tracking Analyses 

Then, we examine the convergence and tracking analyses 
of the proposed CTVFF mechanism with the RLS receiver 
The steady-state MSE between the desired and the estimated 



Fig. 9. BER performance versus the (cycles/symbol) in multipath 
fading channels. We use SNR= 15dB, K= 5, 1500 symbols are 
transmitted and 250 symbols in TR. 



N=15, K=6, f,T=0.0001 



N=1 5. SNR=1 5 dB, f T=0.0001 




6 8 10 12 
(b) number of users (K) 



Fig. 10. SINR performance versus (a) SNR and (b) the number of users 
(K) in multipath time varying channels, f^T = 0.0001. 



symbol obtained through simulation is compared with the 
steady-state MSE computed via the expressions derived in 
Section V. We verify the analytical results (l33T l and (l46t to 
predict the steady-state MSE in time-invariant channels and the 
analytical results ( |40l i and ( |46] l to predict the steady-state MSE 
in time-varying fading channels. In this simulation, we assume 
that four users having the same power level operate in the 
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system. The channel has a power profile given by po = OdB, 
pi = -MB and p2 = -IQdB. We employ 8i = 0.99, 62 = 
0.35 X 10^^ and ^3 = 0.995 for the case of invariant channels, 
and employ 5i = 0.99, 82 = 0.4 x lO'^ and 0.99 for the 
case of time-varying fading channels, where JdT ~ Ix 10^^. 
250 symbols are used in TR. The quantity of £'[q(i)q^(i)] is 
estimated by using the average over 10000 independent exper- 
iments, namely, ( E^'i where = 10000 
and q(i) = wo(i) — wo(i — 1). By comparing the curves 
in Fig. [TT] it can be seen that as the number of received 
symbols increases and the simulated MSB values converge to 
the analytical results, showing the usefulness of our analyses 
and assumptions. 



N=15, K=4 



predicted (time-varying channel) 
-^predicted (invariant channel) 
i^' simulated (time-varying channel) 
simulated (invariant channel) 



N=15, K=4, SNR=15dB 



10' 



-simulatei 
predictec 



LU 
05 



10" 



invariant channel 



500 1000 
number of received symbols 



Fig. 11. Analytical MSB versus simulated performance for conver- 
gence and tracking analyses of the proposed CTVFF scheme. 250 
symbols in TR. 

Fig. [12] shows the effect that the desired user's SNR has on 
the MSB, and a comparison between the steady-state analysis 
and simulation results. We assume that four users operate 
with the same power level in the system. We can see that 
the simulation and analysis results agree well with each other 
1500 symbols are transmitted and 250 symbols are used in TR. 
For the time-varying fading channel, we use fdT = 1 x 10~^. 

VII. Conclusion 

In this paper, we proposed a low-complexity variable forget- 
ting factor mechanism for estimating the parameters of CDMA 
multiuser detector that operate with RLS algorithms. We com- 
pared the computational complexity of the new algorithm with 
the existing gradient-based methods and further investigated 
the convergence and tracking analyses of the proposed CTVFF 
scheme. We also derived expressions to predict the steady- 
state MSB of the adaptive RLS algorithm with the CTVFF 
mechanism. The simulation results verify the analytical results 




fa din 




and 



SNR(dB) 



12. Analytical steady-state MSB versus simulated performance 
the proposed CTVFF scheme in time-invariant and time-varying 
ing channels. 1500 symbols are transmitted, 250 symbols in TR. 



show that the proposed scheme significantly outperforms 
isting algorithms and supports systems with higher loads, 
'e remark that our proposed algorithm also can be extended 
take into account other types of applications. 



Appendix 
^oProo/o/iO 

By multiplying r(i) on both sides of ( fT2] i we have 

-R-\i)Y{{) = A-l(^)R-l(^ - l)Y{i) - A-^(^)k(^)r^(^)R-^(^ - 1) 

(48) 

Rewrite (flOl l as 

k(^) = X-Hi)Kl\i - l)r(i) - A"l(^)k(^)r^(^)R-l(^ - l)r(i). 



(49) 



Based on (081) and ( l49l l we have ( fT9l l. 



B. Derivation of ( 1201 ) 
We note that 

R(i) = A(i)R(i-l)+r(i)r^(i) 

= v{i)v"(i) + \{i)Y{i - l)Y"{i - 1) + - l)Y{i - 2)r^ 

+ • ■ • + \{i)X{i -!)•■• A(2)r(l)r^(l) + - 1) • ■ • A( 

(50) 

For large i, based on (fTST l. ( fTTI l and (fTsT l. we assume that there 
exists a number Ni > 0, when i > Ni, for which we have that 
the forgetting factor A(i) varies slowly around its mean value 

and E[X{N,)] w E[X{N, + 1)] w ■ • ■ « E[X{i)] w E[X{oo)]. 
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By taking the expectation of ( |50l ) we obtain 

E[R.{i)] w E[l + X{i) + X{i)X{i - 1) + . . . + X{i)X{i - 1) 

+ E[X{i)Xit - 1) • . • XiN,)]E[X{N, - l)r{N, - 2) 
+ X{N, - 1)X{N, - 2) • ■ • A(2)A(l)r(0)r^(0)] 
« (1 + E[X{i)] + E^[X{i)] + ■■■ + E'-^'+^[X{i)])E 

(51) 

where the matrix A = £'[A(iV,-l)r(iV,-2)r^(iV,-2)H h 

X{Ni - 1)X{N, - 2)---A(2)A(l)r(0)r^(0)] is determinate. 
Note that < i?[A(z)] < 1, when z — > cx), we obtain 



mm 



1 - E[X{oo)] 



E[r{^)r"{^)]. 



(52) 



Note that, when i — > co, R(i) and R~^(i) converge. 
Thus, we assume linii^oo R"""^!*) ~ ^inii^oo E[Il^^{i)] « 
limi^oo -E"^ [R(i)]. Finally, we obtain ( |20l i. 
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